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1 INTRODUCTION 

1.1 Motivation for Spin Studies 

In recent years there has been a significant interest in the study of spin phenomena. It is 
now widely accepted that the physics of spin phenomena in particle interactions provides 
vital information on the most profound properties of particles: their wave functions, the 
short and long distance dynamics of the quark and gluon interactions, and the mechanisms 
of chiral symmetry breaking and confinement. These systematic spin studies have several 
well defined goals; some of these include: 

• The study of the spin structure of the nucleon, i.e., how the proton's spin state can 
be obtained from a superposition of Fock states with different numbers of constituents 
with spin. 

• How the dynamics of constituent interactions depend on spin. 

• Study the overall nucleon structure and long range dynamics. 

• Understand chiral symmetry breaking and helicity non-conservation on the hadronic 
level. 

The first two goals are related to the constituent short-interaction dynamics, while the last 
two concentrate on hadronic long-range (non-perturbative) dynamics. Experiments designed 
to investigate all of these important properties of nucleons can be performed with polarized 
proton beams at many present accelerators. These experiments can also be used to test some 
fundamental assumptions of QCD regarding spin. 

The naive quark model has been successful in predicting most of the gross properties 
of hadrons, such as charge, parity, isospin and symmetry properties in relation to each other. 
Some of the dynamics of particle interactions can qualitatively be understood in terms of 
this model, as well. However, it has been noted, within the past decade, that it falls short in 
explaining the spin properties of hadrons in terms of their constituents. Perturbative QCD 
(PQCD) has also been successful in predicting asymptotic properties of hadronic dynamics 
in the limit of short range interactions. Once again, some of the PQCD predictions regard- 
ing spin have disagreed with data, even in those kinematic regions where PQCD is thought 
to be valid. The data include measurements of analyzing power, which yielded non-zero 
values, in contrast to PQCD predictions and oscillations in exclusive p — p scattering cross 
sections. These are related to the assumption of helicity conservation on the hadronic level. 



Hyperon production experiments from inclusive p — p scattering have shown that polarized 
hyperons can be produced from unpolarized protons. This is totally unexpected from PQCD 
arguments and the assumption of the pointlike constituent dynamics. These phenomena are 
likely related to longer range dynamics and the coordinated rotation of constituents in terms 
of angular momentum. 

Thus, in probing hadronic structure with regard to spin, a number of interesting 
events have arisen, which require modification of the models created to explain hadronic 
matter. Outside of new particle discoveries, spin phenomena have provided the high energy 
community with a great number of surprises, which have forced us to refine our picture of 
elementary matter. 

In this review, I will discuss a particular aspect of this exciting field, namely the 
quest for understanding nucleon spin with regard to the quark and gluon constituents and 
their overall collective motion. Naturally, there is some overlap in the understanding of the 
constituent spin and the overall nucleon spin, since the constituents must be confined to the 
nucleon "bag". Also, since the collective motion of the constituents plays and important role 
in structure and the dynamics of the interactions, some of the analyses will include the non- 
perturbative effects which are characteristic of the long range dynamics. The main thrust of 
this paper, however, is to describe the theoretical and experimental aspects of understanding 
the constituent structure of nucleon spin in fair detail. 



1.2 The Spin Crisis 

Consituent spin structure studies involve inclusive processes such as the deep inelastic scat- 
tering (DIS) of leptons on nucleons. There are two crucial reasons why deep-inelastic scat- 
tering is the primary tool for probing hadronic structure: (1) the pointlike nature of the 
leptons allows the probing of hadrons so that constituent interactions can be ignored, and 
(2) the theoretically calculable quantities which characterize the constituent structure can 
be directly related to the measurable cross sections. In other words, the cross sections 
are factorizable into a calculable part and a measured part. Historically, the deep-inelastic 
scattering of electrons first proved the existence of point-like constituents in the nucleon. 
The first spin structure measurements of the proton were performed in the 1970's when the 



SLAG- Yale collaboration measured the spin structure function gi{x). |]1| In the late 1980's, 
the European Muon Collaboration at CERN extended the kinematic range of the DIS 
measurements using a polarized muon beam. The extrapolation of the polarized structure 
function Qi to lower x (parton fractional momentum) led to implications that, although the 
Bjorken sum rule of QCD was satisfied, the Ellis- Jaffe sum rule, based on a simple 
quark model with unpolarized strange quarks, was violated. This created a controversy, 
labeled the "spin crisis" which heightened the interest in spin phenomena. Since then, a 
flurry of theoretical and experimental work has been performed to address this "crisis" and 
further investigate the spin properties of the lighter hadrons. 

Initial analyses of the EMC data were motivated by comparison to the naive quark 
model and the Ellis- Jaffe sum rule. This sum rule was based on the assumption of an 
unpolarized strange sea, motivated by unitarity arguments, and will be discussed in detail 
later. The "spin-crisis" arose since the data disagreed considerably with the naive quark 
model, which predicted that all of the nucleon spin should be carried by the valence quarks, 
and the Ellis- Jaffe sum rule. The EMC data led to the implications that: 

• (1) the total quark content of the proton spin is small and possibly consistent with 
zero, 

• (2) the predictions of the E-J sum rule are severely violated. 

There were a number of proposed "fixes" for these problems. Among them were: 

• (1) introduction of the gluon anomaly, which can modify the spin content of each quark 
flavor, so that the measured and calculated values differ by a calculable factor, 

• (2) proposing that the spin carried by the polarized gluons was extremely large ((AG) ~ 
6), 

• (3) assuming that the corresponding orbital motion of the constituents was very large, 

• (4) assuming that the polarized sea was large and negative, canceling most of the 
valence quark contribution to the spin, 

• (5) abandonment of the quark model in favor of alternatives, such as the Skyrme model. 



All of these were motivated by the prevailing notion that the naive quark model with 
an unpolarized strange sea was an accurate model of the nucleons' spins. Analyses of the 
more recent data (from 1992 on) are more inclined to accept the notion that the naive quark 
model must be modified to account for spin contributions of the constituents and that the 
sea is polarized negatively with respect to the valence quarks. This virtually eliminates the 
necessity of making the polarized gluons or the angular momentum extremely large. There 
may be other considerations which affect these quantities in the opposite way, such as nega- 
tively polarized gluons and the possibility of rapidly growing structure functions at small x. 
Thus, in response to the recent DIS data, previling theories have reduced the possible range 
of values of the constituents' spin contributions and have made us reconsider the important 
elements which modify these contributions. 

Recently, much attention of the HEP community has focused on the physics of polar- 
ized hadron interactions. In the last couple of years, experimental groups at SLAC, CERN 
and DESY have improved statistics and lowered the systematic errors in these DIS exper- 
iments. [P, ^, H This data has allowed us to draw several conclusions about the nucleon 
spin structure in the kinematic region of lepton momentum transfer, Q^ — 2 ^ 10 GeV^ 
and Bjorken x ~ 0.002 -^ 0.7. The analysis of these data in the framework of perturba- 
tive QCD provides information on longitudinally polarized parton distributions Agj(x,(5^), 
interpreted as the differences of probabilities q~^^~{x,Q^) for finding partons of the type i 
with spin parallel/antiparallel to the spin of the parent nucleon. The results have required 
modification of the earlier quark models of nucleon spin and have created other controversies 
regarding the amount of spin carried by the strange sea quarks and the gluons. [§] Thus, 
although there has been significant progress in understanding nucleon spin structure, present 
data have left unanswered questions regarding: (1) flavor dependent spin structure, (2) re- 
lations between the polarized and unpolarized quark and gluon structure functions and (3) 
choice of factorization of the quark distributions and the role of the gluon anomaly in their 
determination. 

Recently, the E704 group at Fermilab has analyzed the data from their pp and pp ex- 
periments for jet production cross sections. [|T^] The results have implied that the polarized 
gluon distribution is limited in size, but has not put a strict value on this limit. In all of the 
inclusive experiments done this far, we have gained some knowledge of the constituent spin 
content, but it is clear that much more work has to be done in future experiments to obtain 




Figure 1: Diagram of Deep-Inelastic Scattering, showing relavent momenta. 

a complete understanding of this important problem in physics. Since the initial studies of 
constituent spin structure have relied heavily upon deep-inelastic scattering experiments, it 
is useful to outline the connection between theoretical calculations and the measured cross 
sections. 



1.3 Deep-Inelastic Scattering 

1.3.1 Formalism 

Since the pointlike lepton is incident on the hadron target at very high energy, it effectively 
strikes one of the quarks within the hadron. From the angular distribution of the scattered 
lepton, the momentum information can be inferred. The polarized deep-inelastic scattering 
process (DIS) can be represented by the diagram in Fig. 1. 

In this figure, the initial proton (or hadron target) and lepton 4— momenta are given 
by p'^ and k'^, respectively. The virtual exchanged photon has momentum q and the outgoing 
struck quark and lepton momenta are p^ and k'^', respectively. The cross sections are usually 
written in terms of the invariants: Q^ = — g^, u = E — E' and x = Q'^ /{2Mv). Here, Q^ 
is the momentum (squared) transfered to the nucleon by the lepton, v is the energy lost 
by the lepton in the collision and thus transferred by the photon to the quark, and x is 
a dimensionless variable, called the Bjorken scaling variable, which is a measure of the 



fraction of momentum carried off by the struck quark. [O] Cross section information and 



the corresponding theoretically generated structure functions are represented in terms of 



either the variables v and Q^ or x and Q^ . 

Let us briefly review the formalism of deep inelastic polarized electron scattering on 
a polarized proton. |T2[ Since the thrust of this review is to study nucleon spin structure via 
the constituents, we will cover only the distribution part of the interaction, which involves 
the structure functions. The final-state fragmentation functions are not considered here, 
but the reader can refer to recent treatments of spin measurements of the fragmentation 



functions done by Anselmino, et. ai, |jT^ and Mulders. [|l4| 



In the distribution function case, the hadronic tensor carries all of the relavent infor- 
mation about the target nucleon. This tensor, labeled Wf^^,, will depend on the momenta 
p and q and on the covariant pseudovector of the proton spin, s. Taking into account 
translation invariance, the hadronic tensor can be rewritten as a Fourier transformation of 
the matrix element of the commutator of two currents representing the dynamics of the 
interaction: 

W,Ap,q,s) = -^J^e^'Hp,s\[MO,M0)]\P.s) = W^^ + zWl^K (1) 

where \p, s) denotes the proton state with momentum p and spin s. The hadronic tensor 
contains all of the appropriate information necessary to extract the spin-averaged (unpo- 
larized) and spin-weighted (polarized) structure functions in inclusive scattering. For the 
purpose of separating the unpolarized and polarized information, we can separate the sym- 
metric [S] and anti-symmetric [A] parts of the hadronic tensor. The functions PFJ^fl and 

ly [^1 are determined by the discontinuity of the symmetric and antisymmetric parts of the 
forward Compton scattering amplitude: WJ^^'^' = ^ImT'f'"^!. Imposing parity (P), charge 
conjugation (C) and current conservation {q^W^y = 0) we can extract the forms: 

W^? = {-9,u - ^-^)W,ix,Q') + ^(p, + ^q,){Pu + ^qu)W,{x,Q') (2) 

and 

^^' = -^,e,yx.q'[M's^G,{u,Q') + {Mp-qs^-s-qpnG2{iy,Q% (3) 

This last relation provides the opportunity to apply the helicity amplitude formalism which 
is especially useful in the case of targets with spins greater than |. For a spin-| target 
there are only four independent helicity amplitudes in forward Compton scattering and 



hence four independent structure functions. The spin-dependent structure functions Gi 
and G2 determine the asymmetries which depend on the initial state lepton and nucleon 
polarizations. 

In the Bjorken limit {v and Q^ -^ 00, x constant) scaling is valid and the functions 
Wi, W2, Gi and G2 should depend only on x, up to the logarithmic corrections: 

MWi{u,Q^)^Fi{x), 
iyW2{iy,Q^)^F2{x), 
M'uG,{u,Q^)^g,{x), 
Mu^G2{iy,Q^) ^ g2{x). (4) 

Studies of deep inelastic scattering processes in the Bjorken scaling limit provide us with the 
knowledge on the structure of the target nucleon. In the Bjorken scaling limit, the symmetric 
and anti-symmetric parts of the hadronic tensor are: 

The functions Fi and F2 are the structure functions extracted from unpolarized DIS, while 
Qi and g2 are extracted from polarized DIS. These reveal the pointlike interaction between 
the hard photon and the constituent partons within the hadron. 

In the operator product expansion of QCD, the local operators are quark and gluon 
operators, characterized by definite values of twist. The hadronic electromagnetic current in 
the quark model has the form: 

uo=j:4'i(OiMo- (6) 

<? 

Here q{C,) denotes the quark field. The matrix element between the free quark states is 

(g|JJ|g)~c,(g|0,|g)+c,(g|0,|g), (7) 

where indices q and g refer to the quark and gluon operators. Since the electromagnetic 
current is a quark operator the left hand side and first term in the right hand side are of 



order (a;s)°, whereas the matrix element {q\Og\q) is of order (tts)^, because there are at least 
two gluons entering the operator Og. Thus, in general one can neglect the gluon operator. 
The most important operators in OPE are the operators with the lowest possible twist. The 
twist-two operators correspond to the vector and axial vector currents for the process of 
the virtual Compton scattering which determine the hadronic tensor W^v These operators 
provide the finite contribution to the structure functions in the deep inelastic limit while 
the contribution of twist-three operators are suppressed by MjQ. Transverse spin structure 
functions, for example, are of the twist-three type. 

The structure function g\{x, Q"^) is related to the longitudinal polarization of the pro- 
ton spin with respect to its momentum, i.e. its helicity. In the parton model, gi effectively 
measures the quark helicity density. The second spin-dependent proton structure function 
g2{x,Q'^) is related to transverse polarization of the nucleon spin. It has been measured at 
SLAG (E143) in a limited x region. The analysis based on the OPE does not depend on the 
type of nucleon polarization (longitudinal or transverse). The analysis of the second struc- 
ture function g2{x,Q^) may be performed similar to the function gi{x,Q'^). The significant 
difference between the two cases is that gi{x, Q"^) receives a contribution only from twist-two 
operators whereas g2{x,Q'^) gets contributions from both twist-two and twist-three opera- 
tors simultaneously. A few theoretical results have been obtained for the function (72 (x), but 
simple partonic interpretation is only possible for the twist-two operator contribution. The 
twist-three contributions are not well understood. The following relation exists between the 



functions gi and the twist-two contribution to (72: H^ 

g2ix,Q')= f'%,(y,Q^)-g,(x,Q'). (8) 

Jx y 

The above relation is used to calculate g2{x, Q^) from gi{x', Q"^) aX x' > x in the framework 
of the parton model with free on-shell partons. However there are no reasons to neglect the 
contributions of twist-three operators at low enough Q^ and therefore g2 can be represented 
as follows: 

g2{x,Q^) = g2{x,Q^f^ + g2{x,Q^f\ (9) 

where the first term on the right hand side of Eqn. (9) is provided by Eqn. (8). The twist- 
three operator contributions (72(3^, <5^)'^^ depend on the effects of quark-gluon interactions 
and quark masses. 



Experimental measurements oi g2(x, Q"^) provide a direct way to study the magnitude 
of twist-3 contributions and accurate data with polarized proton beams would be crucial 
in the resolving of the question of the magnitude of these contributions. Data could also 



provide a test for the Burkhardt-Cottingham sum rule: [16 







^2(x) = 0. (10) 



The success or failure of Eq. 1.10 depends on its long-range behavior which, in the simple 
models, is such that this sum rule is satisfied in perturbative QCD. Thus, valuable informa- 
tion regarding both perturbative and non-perturbative processes can be obtained by testing 



this sum rule. Recent measurements of g2 by the SMC and E143 [T^ experiments have 
verified the condition that A2 < \/R {R = ci/aT) so that the leading twist terms dominate 
and that the Burkhardt-Cottingham sum rule is valid to within experimental errors. There 
are a number of recent treatments of the role of g2 in DIS with regard to transverse spin 



content and chiral symmetry breaking. [18] This topic will not be discussed here 



The above discussion has shown how the hadronic tensor can be written in terms of 
the spin-averaged (unpolarized) structure functions Fi and F2, and the spin- weighted (polar- 
ized) structure functions gi, g2- It is these quantities which can be extracted from the cross 
sections in various DIS processes. The phenomenology of DIS starts with combining the 
incident lepton tensor with the hadronic tensor and relating this product to the measured 
cross sections. 



1.3.2 Phenomenology 

The polarized lepton tensor L^^, being pointlike, can be written explicitly as: 



^fiu — ^ 



'^fi'^u + '^u'^u 9fJ,iy'^ ' "^ ^^uj "I" '^^fiupo-Q S , (11) 



where the antisymmetric last term is absent in the unpolarized case. The product of L^i^ 
and W/j,^ is directly related to the cross sections measured in the process of DIS. 

Since cross sections are the most fundamental measured property of elementary par- 
ticles, all of the theoretical quantities characterizing spin must be extracted from some 



combination of the polarization cross sections. [0 The following differential cross section 
relates to the lepton and hadron tensor product as: 






(12) 



This connects the structure functions in the Bjorken limit directly to the measured cross 
sections, both in the unpolarized and polarized cases. These differential cross sections can 
be written as in terms of any of the kinematic variables. For example, if the energies and 
angular distributions of the leptons are measured, the unpolarized and polarized structure 
functions are related to the spin averaged and spin weighted differential cross sections by: 



d^ai 



dVtdE' 



+ 



t/V( 



2 771/2 



dVtdE' 



MQ4 



2sin2(^/2)Fi(x,Q2) + (MVz^)cos2(e/2)F2(x,Q2) , (13) 



and 



rfV( 



d^ai 



Aa'E' r 



dVtdE' dVtdE' Q'^Eu 

where dVt is the angular distribution of the leptons 



{E + E' cosie))gi{x, Q^) - {2xM)g2{x, Q' 



(14) 



In nature, asymmetric behavior often occurs between the interactions of particles 
whose spins are aligned and interactions of particles whose spins are anti-aligned. Thus, it 
is convenient to define an asymmetry as: 



A 



a(++) + (T(+- 



where the ++ refers to the aligned spins and H — to the anti-aligned spins, without regard to 
the direction of polarization relative to the momentum of the interaction. There are distinct 
advantages to defining the asymmetries in this way. Theoretically, A is a ratio of cross 
sections, so all of the normalizations used to calculate cross sections cancel, leaving only the 
key parameters of the interaction. Experimentally, the difference in the numerator tends to 
cancel some of the systematic errors, making A a fundamentally more sensitive measure of 
the physical variables. 

There are asymmetries for both the longitudinal and transverse spin relative polar- 
izations. The lepton-hadron asymmetries are given in terms of the cross-section as: 



A||(x,Q' 



o- 
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(15) 



where the arrows refer to the relative longitudinal spin directions of the beam and target, 
respectively. 

The structure functions can be related to the absorptive cross sections, ai and as for 

2 2 

virtual photons with helicity projections | and |, respectively. The asymmetry components 
Ai and A2 can be written in terms of the corresponding cross sections along with the inter- 
ference between the longitudinal and transverse polarizations, aj. The parallel asymmetry 
in terms of these components can be written as 



A\\=DiAi + riA2), 
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ai 
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(16) 



The kinematic factors are: D = i,^^ and r] = 3]^_^_ \. , where y is the fraction of energy 

lost be the lepton in the lab frame and R = a 1/(^X1 the ratio of transverse to longitudinal 
photon-nucleon (7 — A^) cross sections. In most DIS experiments, the kinematics are such 
that D is known for each experiment and t] ~ eQ/ s is very small. The asymmetries can be 
linked to the structure functions gi and (72 by the following sequence of relations. First, the 
longitudinal asymmetry is factored into a combination of the longitudinal j — N asymmetry, 
Ai{x, Q"^) and the transverse 7 — A^ asymmetry, A2{x, Q"^) in the following way: 



A\\ =D{Ai+riA2) 



D 



[1 + 7% 2 ^^ 7(1 

(^1 - 7 92)+ V 



r 



Fi 



Fi 



-{91 + 92] 



(17) 



where 7^ = {2MxY /Q"^. From recent measurements, it has been shown that A2 is small 
[E143, 1995 and 7 is also negligible for the SMC experiments. In addition, it is also assumed 
that the transverse asymmetry A2 is small, since it is bounded by ^JR. As a result, the 
structure function ^ff is neglected in both factors of A\\. Experimental measurements appear 
to imply that Ai is relatively independent of Q^, an assumption which has been substantiated 
by phenomenological analysis of the data. [|l^ However, Gliick, et. ai, [^ point out that 
this may not strictly be a valid approximation at low-x and moderate Q^. Assuming that this 
is a reasonable approximation, then, in the Bjorken limit, Fi = 2xF2, and thus, Ai ^ Ah/D 



with 



(/i(a;,Q2)^Fi(x,Q')-A(a;) 



2x(l + R) 



Information about the polarized quark distributions can be extracted directly from this 
asymmetry by 



Eietqi{x) 
where Cj is the charge of each quark flavor and for each flavor, i: 



(19) 



(20) 



Using the OPE it can be shown that the first moment of the proton structure function 
gi{x, Q"^) is determined by the following matrix elements of the axial- vector current: 



dx gi{x, Q^ 



^AuiQ') + lAdiQ') + Ias{Q') 



1- 



asiQ' 



Oial 



O 
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(21) 



where 



Ag(/i^)s^ = (p,s|(g7i.75g)L2b,s). 



(22) 



The term /i^ is the relevant mass scale or the renormalization point for the axial-vector cur- 
rent operator. The functions Aqi^Q"^) are related to Ag(/i^) by the QCD evolution equations 
that will be discussed later. Note that another leading twist operator, the vector current 
qit^q, provides the finite contribution to the symmetrical part of the hadronic tensor IV^f'. 

Combining the form of Eqn. 21 for the proton and neutron, we arrive at the polarized 
version of the Bjorken sum rule, which relates the first moment of the difference between the 
proton and neutron structure functions, g^ — g\ to nucleon beta decay: 



1 \ r\ 1 as(<5^) 

{g{{x) - gl{x)\ dx = - [Autotai{x) - Adtotai{x)] dx = -5'a(1 h h.o.c), (23) 

o JO o vr 



where gA is measured in nucleon beta decay and h.o.c. refer to calculated higher order QCD 
corrections. This and other sum rules will be discussed later. 



The higher order corrections to DIS structure functions have been analyzed by various 
groups. In particular, the higher twist corrections to the proton and neutron functions by 
Stein, et. ai, pl| and were found to be Q^ dependent, but very small, even at the lower 



Q^ values of the data. The higher twist corrections to the BSR were consistent with zero. 
Meyer-Hermann, et. ai, calculated the twist-4 contributions to gi using the IR-renormalon 



method |22| and obtained a result of ±0.017 GeV^/Q^- (See also Boer and Tangerman p3| ) 
The twist-2 and twist-3 contributions to the BSR are also suppressed by a factor of l/Q^- 
The higher order QCD corrections will be discussed with the sum rules in section III, since 
they play a major role in the phenomenology of polarized DIS. 

The rest of this review is structured as follows. In section II, we discuss theoretical 
models for the spin-averaged (unpolarized) distributions and the spin-weighted (polarized) 
distributions. The important aspects of how the theoretical models can be compared to data 
is discussed in detail throughout the section. The experimental data and comparison of the- 
ory and experiment is covered in section III. Physical consequences of this comparison and 
the open physics questions are addressed. Finally, we outline a set of experiments which can 
be used to address these questions and further refine our understanding of the constituents' 
contributions to the spin of nucleons. 



2 THEORETICAL BACKGROUND 

2.1 Unpolarized Distributions 

2.1.1 Formalism 

In the last section, we discussed the unpolarized and polarized structure functions for par- 
ton distributions in the context of the operator product expansion and their extraction from 
deep-inelastic scattering (DIS) data. In this section, we will outline the details of modeling 
the constituent distributions and discuss the possible relations between the unpolarized and 
polarized versions. The important aspects of factorization and evolution will be covered, as 
they are crucial to relating the theory to experimental data. 



Due to the models of the polarization mechanism, |24, 25] and the possible relation 



of the polarized to the unpolarized distributions, knowledge of the unpolarized distributions 
is useful before one analyzes the polarized case. Naturally one can extract the x-dependent 
polarized distributions from the data directly, |^U|, |2^ but the former method has a more 
direct theoretical connection. Normally, one assumes an SU(6) symmetric wave function 
for the proton, from which the polarized valence distributions can be found from a suitable 
parametrization of the unpolarized ones, Uv{x) and dv{x). Most early models of the polarized 
sea and gluon distributions also used this assumption as a starting point. ||2^, |29l 



For purposes of studying the quark properties of nucleons, it is convenient to dis- 
tinguish between the "valence" and the "sea" quarks of the proton. The valence quarks 
carry the main quantum numbers of the nucleon. For the up and down quark densities in a 
nucleon, we can write 

u{x, Q^) = Uy{x, Q^) + Us{x, Q^) 

dix,Q^)=d,ix,Q^) + ds{x,Q^). (24) 

The valence distributions are the flavor nonsinglet components of the quark distributions in 
the proton and are normalized so that 



/ dxUy{x, Q"^) = 2 

/ dxd^{x,Q^) = 1. (25) 



The neutron, of course, exchanges the roles of the up and down quarks, since it consists of 
the "udd" combination of valence quarks. All other flavor components will be included in 
the sea. In general, these distributions will be denoted by their flavor quantum numbers; 
u{x,Q^), d{x,Q^), Us{x,Q^), ds{x,Q^), s{x,Q^), s{x,Q^), etc. 

For unpolarized distributions, there presently exists a large body of phenomenological 
knowledge concerning the various parton densities, which have been compiled into various 
sets of models, which will be discussed shortly. The differences between the different param- 
eterizations of the unpolarized distributions are usually minor compared to the uncertainties 
in data from which they are generated. 



2.1.2 Evolution 



The Q^ evolution of the unpolarized (spin-weighted) structure functions is governed by the 
GLAP equations, [BO] which in the leading logarithm approximation take the form 



— [Qv{x, t)] = Pgq{x) O Qy{x, t) 

- [Q(x, t)] = Pqg{x) ® Q{x, t) + 2NfPqG{x) ® G{x, t) 



dt 



where, 



G{x, t)\ = Pcgix) ® Q(x, t) + Pgg{x) ® G{x, t), 
Q{x,t) =^a;gi(x,t). 



ll-fiV, 



In 



2 / a2n 



and 



ln(Q7A^)/ln(Q^/A 



Pix)0Qix,t)= [ ^P(z)Q(-;t). 
Jo z \z / 

The probability kernels, Pij{x), for three quark flavors are given by 



^qq{X) 



1 — X 



+ 26{l-x) 



Pqcix) = -[x' + {1 



X] 



Pg,{x) 



Pgg{x) = 3 



l + (l-x)^ 



X 



X 



1 — X 



h x(l -x) + -5(1 - x) 

X 4 



(26) 



(27) 



f28) 



(29) 
(30) 
(31) 



The (1 — x)+ distribution renormalizes the kernels to avoid divergence at x = 1 and is defined 
by: 



h 71 



dxj{x) ^ ri^J(x)-f{l) 



(1-x). 



1 — X 



(32) 



These equations imply that even if there is no initial gluon distribution generated at 
the nonperturbative level, a positive G{x) will be generated by quark Bremsstrahlung, since 
Paqix) is everywhere positive. 

The next-to-leading order (NLO) splitting functions have been calculated for the un- 



polarized distributions. |31| The separation of the non-singlet (related to the valence distri- 
bution) and the singlet evolution provides a convenient motivation for the separation of the 
valence and sea quarks in the analysis. Thus, the evolution equations allow us to determine 
the Q^ dependence of the valence and sea densities. We can also investigate the small-x 
behavior of the valence quark and gluon evolution, which play a significant role in models 

of the nucleon. There seem to be differences in some of the NLO evolution codes used to 
generate unpolarized structure functions at small- a; , primarily due to truncation errors in the 

NLO terms, i.e. in the NNLO terms. [^ These manifest themselves most dramatically at 

low X and high Q^. This is a minor difficulty relative to other uncertainties of the structure 

functions at small-x, as will be discussed later. 



2.1.3 Existing Distributions Prom Data 

Unpolarized parton distributions are determined from global analyses of data for a large 
range of processes in as wide of a kinematic range as possible. There have been various 



groups who have generated these distributions. ||3^, |3^, |35| The improvements in these dis- 
tributions have come from more precise data in a wider kinematic domain. [^] Theoretical 
progress in calculating higher order contributions, both logarithmic (to a^) and power law in 
Q^ (higher twist), have increased the accuracy to which these distributions can fit the data. 
These also help obtain a more rigorous test of various aspects of QCD. The differences in 
the unpolarized distributions are as follows. 

The "Martin, Roberts and Stirling" (MRS) models P3 extract the unpolarized distri- 



butions from data, normalized to a fixed Qq value, normally around 4.0 GeV^. The valence 
and each fiavor of sea quarks are separately parametrized and the excess of anti-down over 
anti-up quarks is built in to the analysis. The parton distributions have the general form: 

xD{x) = Nx"" {1-xf {1 + ax^ +bx), (33) 



with the unkonown factors determined by various data. The gluon density is chosen to be 
finite at X = and the sea has a Regge type behavior at this hmit. 



The "Gliick, Reya and Vogt" (GRV) distributions [Q derive x and Q^ dependent 



distributions from the data. The valence parametrization is done separate from the sea and 
the up and down fiavors along with the gluons include a logarithmic dependence at small-a;. 
The strange sea is modeled separately from the lighter fiavors. The form for the valence 
distributions is: 

xQ{x, Q2) = Ar(Q2)3,a(Q2) ^i _ ^y{Q^) [1 + a(Q2)^i ^ b{Q^)x] (34) 

and the general form for the other unpolarized distributions can be written as: 
xSix, Q2) = iV(Q2)a;"(Q^) (i _ j^fiQ') [i + a{Q^)x + b{Q')x'] [in(-) + /(In(-), Q^] . (35) 

Most of the parameters are Q^ dependent, evolved both in leading and next-to leading orders. 

The CTEQ distributions B^, like the MRS, are normalized to a fixed Qq. They have 
the general form: 

xQ{x) = Aox^' (1 - a;)^2 (^ ^ Agx^^), (36) 

with the parameters determined by fits to data. The MRS and CTEQ distributions are 
evolved from the Qq value by using the NLO GLAP equations. [^, ^ All of these dis- 



tributions are useful for calculating observables characteristic of the unpolarized scattering 
experiments. Further, if suitable assumptions are made about the polarized distributions, 
then their x dependence can also be extracted similar to these parametrizations. 



2.1.4 Unpolarized Distributions at Small-x 

Differences in the above distributions for the most part lie in the small-x region. Recent data 
from HERA has shed light on the small-x behavior of the parton distributions, but questions 



of the parametrization and evolution of this behavior remain controversial. |^ Since the 
experimental errors in most polarized DIS experiments are the largest in this region, there 
is no apparent advantage that one unpolarized distribution has over the others in terms of 
generating the corresponding polarized distributions. Most of these distributions build in a 
parametrization for the difference in anti-up and anti-down quarks, first discovered by the 



New Muon Collaboration IBSi measurement of the Gottfried sum rule at CERN ||39ll and 



later confirmed by the NA51 measurement of the pp/nn asymmetry in Drell-Yan. |^ For 



further discussion on the asymptotic behavior of unpolarized distributions at small-x, see 



Webber pT| and Lopez, et. ai, H3] The polarized distributions' small-x behavior will be 



discussed in the section covering the polarized sea. 

2.2 Polarized Parton Distributions 

2.2.1 Factorization 

The usual application of the QCD-parton model assumes that the appropriate distributions 
are measured in one set of processes. Predictions are then made for other processes by 
using factorization. For the purpose of estimating the impact of possible experiments, it is 
convenient to have model distributions which represent an informed guess concerning the 
allowable range of distributions. 

The two basic ideas which will help us in constructing these models are: 

1. There is a strong connection between the spin of the valence quarks and the spin of 
the proton. 

2. There are spin-dependent forces between the constituents in the proton which influence 
the shape of the spin-weighted distributions. 

A successful description of the baryons within the framework of the "naive" constituent 
quark model suggests strongly that, in some approximation, a large part of the spin of the 
proton is associated with its valence quarks. 

When extracting information from hard-scattering processes, such as DIS, it is nec- 
essary to write the measured cross sections as a product of factors calculable in PQCD 
and non-perturbative (usually extractable from data) terms. Such a process is known as 



factorization. Factorization theorems exist for many processes including: [43 



Inclusive DIS 

Semi-inclusive production of jets, heavy quarks 

Drell-Yan: lepton pair production from hadron-hadron interactions 



• High pt processes: jet production; hadron production; direct photons 

• Hadron-hadron to heavy quark inclusive production 

• e"'"e~ to jets 

• Elastic processes 

The importance of factorization is that it separates the short distance behavior of the 
cross sections, calculable in perturbative QCD and the long range non-perturbative physics. 
The non-perturbative part should, in principle, be guage- independent, to ensure that it is 
measurable. In principle, one can define factorized quantities which are gauge dependent, 
but it is desirable to have them gauge- independent, since the uncalculable quantities must 
be measurable. In polarized DIS, there are two schemes which are used to define the quark 
spin densities: (1) the gauge-invariant and (2) the chiral-invariant schemes. Since chirality 
is involved in the spin-dependent processes, neither of these is inherently superior over the 
other. 1^ A detailed discussion can be found in Cheng's review. |^ 



Gauge- invariant scheme: EB 



According to an operator product expansion analysis (OPE), (see section I) there is 
no gluonic operator contributing in leading order to Qi. The hard-gluonic contributions and 
the quark spin densities are k± factorization dependent. The factorization scheme which 
respects the OPE is such that the hard gluons do not contribute to the quark spin densities, 
but the hard gluon component of spin does perturbatively generate a negative sea polariza- 
tion due to the axial anomaly in the triangle graph for j^ between external gluons. Thus, 
the anomaly term becomes part of the photon-gluon cross section and we would expect to 
see a correspondingly large negatively polarized sea, due to the gluon Bremsstrahlung which 
creates it. Chiral symmetry is broken in this scheme. 



Chiral-invariant scheme: ||47I 



In this scheme, the polarized quark densities are factorized into a gauge dependent 
piece and the gluon axial anomaly is explicitly separated out. Here, the measured quantities 
extracted from data and sum rules are modified by the anomaly term. This accounts for the 
observed differences between the first moments (integrals) of Fi and gi. The gluon anomaly 
contributes directly to the quark spin densities via the term T^Q"^) added to each of the 



quark flavors. Large negativity of tlie sea is not required, but a smaller polarized sea could 
suggest a larger polarized glue for the naive parton model for the Ellis- Jaffe sum rule to be 
valid. Recent data indicate that this is likely not the case. 

The key differences in the two schemes lie in the k± factorization of the quark spin- 
density and the hard photon-gluon cross section. Then, the factorization prescription is 
determined by choice of the ultraviolet cutoff. In any case, it has been shown that gi is 
independent of which scheme is used. In both schemes, the polarized gluon density plays a 
role in the contribution of each quark flavor to nucleon spin. Should AG be small for the Q^ 
values of present data, the flavor dependence of spin contributions will be approximately the 
same in both schemes. Recent treatments by Goshtasbpour and Ramsey |^, ^ extract the 
flavor dependent information from data using both factorization schemes by setting AG = 
in one of the models considered. These results will be outlined in a later section. Thus, the 
controversy remains as to which of these is more appropriate in explaining the constituent 
spin problem. Further experimentation will determine whether it is the sea or glue that is 
large enough to explain why the quark contribution to the nucleon spin is smaller than unity. 

The choice of a factorization scale, the scale at which the distributions will be used 
in the calculation of various inclusive processes, [^ is consistent in all of the experimental 
and theoretical analyses. 



2.2.2 Sum Rules 

After factorization, the next major step in extracting spin parton densities from the data is 
by use of the sum rules. In section I, we discussed the hadronic tensor and its relation to 
the cross sections of DIS. The targets in DIS are generally characterized by sets of conserved 
quantum numbers. These can often be built into the analysis by forming combinations of 
the cross sections which characterize these conserved quantities. In this way, theory and 
experiment are directly compared. These combinations are referred to as sum rules. Those 
sum rules which involve the strong coupling naturally will have higher order logarithmic 
(perturbative) corrections and since most are valid over a wide kinematic range will also 
have higher twist [Q"^ dependent) corrections. In DIS, the flavor dependent parton spin 
information can be extracted from the data, provided certain sum rules are assumed to be 



valid. There have been recent reviews on the parton model sum rules, [Q but those which 
are most pertinent to the polarized distributions will be discussed here. 

For a proton with its spin aligned in the +z direction (along its momentum), we can 

impose the parity invariance of the strong interactions and define the spin-weighted densities 
as 

= V,{x,Q')q\x,Q'), (37) 

where, for a given fiavor of quark, the factor rii{x, Q^) is called the polarization of that fiavor. 
For an arbitrary direction of spin, the component of spin along the z (momentum) direction 
is the helicity. The spin-averaged (unpolarized) distributions are just the sums of these 
helicity states. If we consider the proton wave function as characterized by momentum p^ 
and spin s^, the polarized distributions integrated over all x can be represented in terms of 
the Dirac matrices 7^^ and 75 by: 

(AqiS^") = {ps I q-f^j^qi \ ps)/2m, (38) 

where m is the mass of the particle. The related axial-vector current operators, A^, are 
members of an SU(3)/ octet, whose non-zero elements provide relations between the polar- 
ized distributions and data. The non-vanishing matrix elements of these operators define 
measurable coefficients, a^, which provide some of the sum rule constraints used to extract 
parton spin information from the data. These are defined by 

{ps I 4 I ps) = s,a' (39) 

where the a'' are non-zero for k =0, 3 and 8. The matrix elements, a^ are determined from 
weak decays of processes where fiavor changes occur. 

The 03 matrix element, measurable in nucleon beta decay, occurs in the Bjorken Sum 
Rule (BSR), which is based on isospin invariance. It is considered to be fundamental to 
QCD and its validity has been the basis for much of the recent QCD analysis on polarized 
DIS. [^, ^, 1^ The polarized version of the BSR can be written as: 

P-r^ f dx {gl - = ^(1 + a—), (40) 
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where af^^ are higher order logarithmic corrections, calculable in QCD. To 0(af) these are: 
^corr ^ (^) ^ 3,5833(^)2 _^ 20.2153(— )3 + 130(— )^ (41) 

TT TT TT TT 

where the last term is estimated. ||5^ Ellis, et. al, have used the method of Pade approx- 
imates to estimate the higher order corrections to the BSR. |Q Their result is consistent 
with the term quoted above. 

The matrix element as is determined by the weak decay constants, F and D, which 
are constrained by hyperon decay experiments. In terms of the flavor dependent polarized 
distributions, this can be written as: 

as= f dx {Au + Au + Ad + Ad- 2 As - 2 As) = 3F - D. (42) 

Jo 

As Lipkin has pointed out, [Q since a suitable hyperon model does not yet exist, one must 
be careful about imposing the SU(3)/ symmetry and this constraint on the polarized sea. 
Some analyses of the recent DIS data do not rely heavily on this constraint, but use it to 
narrow the relative size of the flavor dependent polarized sea distributions. |^ 

The current Ag is determined by hyperon decay and its eigenvalue ag is related to the 
polarized distributions by: ag = {[Autotai + Adtotai - 2Astotai]) ~ 0.58 ± 0.02. Finally, Oq is 
related to the total spin carried by the quarks in the proton. It provides one part of the 
spin- 1 sum rule, known as the Jz = ^ sum rule. Thus, in terms of the parton distributions. 



ao 



/ dx (Am + Au + Ad + Ad + As + As) = Aqtot- (43) 
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The approximation above is due to the gluon anomaly, to be discussed later. Then, the axial 
currents are related to the structure function Qi in the anomaly-independent form: 

ao = 9(1 - a^T' f\ii^) dx - ]as - ^a^ ^ (Ag^i). (44) 
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The polarized distributions are related to the orbital angular momentum of the constituents 
by the Jz = | sum rule. If we define AG = G+ — G_ and the total orbital angular momentum 
of the nucleon constituents about the z— axis as Lz, then the Jz = \ sum rule is: 

Jz = -Aqtot + AG + Lz = -. (45) 



This represents the decomposition of the constituent spins along with their relative angular 
momentum, L^. 

The Ellis- Jaffe sum rule [Q is less fundamental than the BSR and is based upon the 
assumption that the strange sea is unpolarized. This assumption was based primarily on the 
OZI sum rule, which assumes unitary symmetry for mesons and baryons. |5^ It implies that 



the singlet and octet contributions to the axial-vector currents are equal, i.e., that qq = cls- 
The OZI sum rule has been shown to be valid for mesons, but data imply that it fails for 
baryons. [^ Ellis and Jaffe applied this rule to polarized structure functions to predict their 
averaged values: 

^' ^ ^ 12^ ^3 {F/D) + iy ^^^ 

The + and — signs refer to the proton and neutron cases, respectively. F and D are the 
weak decay constants discussed above. The Ellis- Jaffe sum rule (EJSR) has been widely 
discussed, since the EMC DIS data were published, 0] which indicated a severe violation 
of this sum rule. As previously mentioned, this created the "spin-crisis" which has led to a 
re-assessment of the early models of proton spin. Violation of the EJSR (or OZI sum rule 
assumption for baryons) reduces to the question of how polarized the strange sea is in the 
nucleons. 

Regarding this question, it seems clear from recent data that the strange sea is po- 
larized opposite to that of the valence quarks, but the size of that polarization is still under 
question. There is a theoretical constraint on the polarization of the strange sea, based on 
the positivity of the probability interpretation of the leading order parton distributions in 
the naive parton model. This implies an upper limit to the polarized strange sea contri- 
bution, which may affect its contribution to the Jz = | sum rule. [^ The essence of the 



positivity constraint is that the spin carried by the strange sea is bounded by its momen- 
tum. There may be non-perturbative contributions to this limit, which may affect its precise 
value. This was proposed by loffe and Karliner |^9[ to explain violations of the OZI sum 



rule for baryons. Gehrmann and Stirling [3C] have also pointed out that this limit is only 



valid to leading order, since at higher order, the parton distributions no longer have a strict 
probabilistic interpretation. Most of the existing DIS data appear to violate the limits of 
this bound to some extent. This issue is still under question. 



A recent analysis by Ma [£T| attempts to resolve the EJSR violation by introducing 
a Wigner rotation term as a scale parameter to reconcile the data with the naive quark 
model. This is an interesting idea, but as yet there is no apparent way to strictly test this 
assumption against other models with further experimentation. This may be an area for 
future exploration. 

Another possible mechanism which has been proposed to explain the violation of the 
EJSR is the instanton contribution to the spin-flip mechanism. This is a non-perturbative 
vacuum fluctuation of the gluon field, which gives rise to an effective screening of the valence 
quark polarization by the sea. [^ Predictions for g^ and the Drell-Yan asymmetry are 
given, so that there exist experimental tests for this model. 

There are other sum rules which apply to the transversely polarized structure function, 
g2- The naive parton model predicts that g2{x, Q^) vanishes everywhere. ||50| The Burkhardt- 
Cottingham sum rule [jl6| is somewhat less stringent, namely that 



/ dxg2(x) = 0. (47) 
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This does not imply that g2 is identically zero, in fact recent data indicate otherwise.! 



The Wandzura-Wilczek sum rule ||T^ relates the structure functions gi and g2- 



9i{x) + g2{x)= —gi{y)- (4J 

Jx V 



ix y 

This was obtained by an OPE anaysis of quarks in the infinite momentum frame. The im- 
plication of this sum rule is that the transverse spin of the proton is carried mostly by the 
sea quarks at small- a; . 

An alternate sum rule for g2 is derived from a field-theoretical framework and is not 
dependent upon the OPE, where the hadronic matrix elements are local operators. This 
"ELT sum rule" [0] relates the valence components of gi^2 and is exact for each flavor. It 
has the form: 

dx X [gX{x) + 2gX{x)] = 0. (49) 





There are experimental tests which would give information about this sum rule, namely, 
semi-inclusive meson production or jet production from unpolarized DIS. 



2.2.3 Q2 Evolution 



The Q^ evolution of the polarized structure functions is analogous to that of the unpo- 
larized functions. The polarized version of the GLAP equations in the leading logarithm 
approximation, takes the form 



- [AQ,(x, t)] = APgg{x) (E) AQ„(x, t) 

J [Ag(x, t)] = APggix) AQ{x, t) + 2NfAPqG{x) ® AG'(x, t) 



(50) 



d_ 
di 



AG{x, t) = APcqix) ® AQ{x, t) + APgg{x) ® AG{x, t), 



where, similar to the spin-averaged (unpolarized) case: 

AQy{x,t) = xAqy{x,t) and 
AQ{x,t) = Y,xAq,{x,t). 



(51) 



(52) 



The variable t and the convolution operation have the same form as in the unpolarized case. 
The probability kernels, APij{x), are given by 



(53) 



AP,,{x) = 


4 
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AP,g{x) - 


1 
= 2 
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APGgix) = -{2 - x) 



(54) 



APgg{x 



\^^] +i3-3x + x^ + x^)-^6{l-x) 
1 — X . 12 



These equations imply that even if there is no initial gluon polarization generated at the non- 
perturbative level, Bremsstrahlung from the valence quarks will generate a positive AG{x), 



since APcqix) is everywhere positive. Furthermore this mechanism will also polarize the 
sea. However, since the first moments of AP^g and APgc vanish in the leading order, 
{AQ{x,t)) = /q dx AQ{x,t) is constant in t (and hence, Q^). 

The NLO splitting functions have been calculated and are found to be renormalization 
scheme dependent. [^ Once a renormalization scheme is chosen, these evolution equations 



determine the Q^ dependence of the valence, sea and gluon spin densities, which effectively 
compare the spin carried by constituents to their momentum. These become crucial in 
comparing the different experimental data and their consequences at different values of Q^. 
Bliimlein and Vogt have carried out the NLO evolution of Qi'^ in the MS scheme. |^6| The 



NLO evolution plays an important role in the small-x behavior of these structure functions. 



371 This will be discussed later. 



The Q^ evolution of g2 has been investigated. [^, [T8| Since these functions are small 
and not directly pertinent to the constituent spins, this will not be discussed in further detail 
here. Instead, we will proceed directly to the various polarized constituent distributions. 



2.2.4 Valence Quark Models 

Fundamentally, we assume that the nucleons are comprised of valence quarks, whose polar- 
ized and integrated distributions are defined by: 

Aq^{x, Q^) = g+(x, Q^) - g^(x, Q^) 

{Aq^{Q^))= f Aq,{x,Q^)dx, (55) 
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where +(— ) indicates the quark spin aligned (anti-aligned) with the nucleon spin. In order 
to construct the polarized quark distributions from the unpolarized ones, we can start with 
a modified 3-quark model based on an SU(6) wave function for the proton. This model is 
based on flavor symmetry of the u- and d-sea and constructs the valence distributions to 
satisfy the Bjorken sum rule. The valence quark distributions can be written in the form: 

2 
Au^{x,Q'^) = cos 6 D[uy{x,Q'^) - -dy{x,Q'^)], 



Ad^{x,Q'^) = --COS 6 Ddv{x,Q^), (56) 

where cos6d is a "spin dilution" factor which vanishes as x — *> and becomes unity as x ^ 1, 



characterizing the valence quark helicity contribution to the proton. pU| , pPI Normally, the 
spin dilution factor is adjusted to satisfy the Bjorken sum rule and to agree with the deep- 
inelastic data at large x. 

Two-body spin-dependent forces have a direct influence on the spin-weighted quark 
and gluon distributions, and with simple assumptions about their parametrizations, Qiu et. 
al. , have derived a form for the valence spin-dilution factor cos Oj:, . This spin dilution factor 
has the form 

cos^b(x) = ^^ ^ [1 + N{Q^) ■ X ■ G(x, g2)]-i (57) 

where it is assumed that xq^{x) « xG^{x) at small values of x. The N{Q'^) factor is a 
normalization term, which is adjusted so that the valence distributions satisfy the BSR. Note 
that at small-x, the spin dilution factor has the form: 



cos 6' 



D 



l + N ■ A„-x 



l-CXc. 



9 



"-i^^-^°'" (5^ 



where Ag and ag are the normalization and small x power coefficients for the unpolarized 
gluon distribution, respectively. 

The integrated polarized structure function, P^"^ = /q gi (x) dx, is related to the 
polarized quark distributions by 

/"^"^ = ^(1 - ari{m)^utotai + m^dtotai + as,,,„,]), 

where the QCD corrections (af^^) are given in Eqn. (41). In terms of the polarized distri- 
butions and the assumptions of a flavor symmetric polarized u and d sea, the BSR can be 
reduced to: 

/ \Au^{x, Q^) - A4(x, Q2)] dx = a3{l- — + h.o.c). (59) 
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Thus, the valence contributions can be determined uniquely by this model. Any of the 
unpolarized distributions in principle can be used to generate the valence quark distributions, 
evolved to the Q^ scales of each experiment. These all agree for x > 0.05, but have subtle 



differences for the smaller x values. The spin dilution factor is determined from the BSR, and 
differences in the unpolarized distributions are compensated by adjusting the normalization 
factor A^ in the spin dilution term. Thus, the valence distributions are not sensitive to the 
unpolarized distributions used to generate them. The consistency of the resulting polarized 
distributions can be checked by comparing them with the value generated for the ratio of 
proton and neutron magnetic moments: 

/ip ^ 2{Au^) - {Ad^) _ _3 

fin 2{Ad,)-{Au,) ^2' ^ ^ 

Using the values (Auy) = 1.00 ±0.01 and {Ad^) = —.26 ±0.01, both the BSR and magnetic 
moment ratio are satisfied. This also yields a spin contribution from the valence quarks 
equal to 0.74 ± 0.02, consistent with other treatments of the spin content of quarks. 0, |7l| 
The quoted errors arise from data errors on qa/ Qv and any small differences remaining in 
the choice of the unpolarized distributions used to generate the Aq^ terms. The original 
analysis by Qiu, et. ai, |^9|] effectively reached the same conclusion. 



By an appropriate assumption regarding the relation between the polarized and unpo- 
larized distributions, an a;-dependent set of polariaed valence distributions can be generated. 
Ma [^ has used a quark-spectator theoretical model to generate an alternate set of polar- 



ized valence distributions. These include Wigner rotation parameters which are fit to the 
data. Thus, the approach for generating the valence terms is different than the above model 
and comparison would be parameter dependent. 



2.3 Models of the Polarized Gluons 

2.3.1 Overview of Models 

The gluons are polarized through Bremsstrahlung from the quarks. The integrated polarized 
gluon distribution is written as 

(AG) = / AG{x,Q^)dx= f [G+{x,Q^)-G-{x,Q^)]dx, (61) 
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where the ±(— ) indicates spin aligned (anti-aligned) with the nucleon, as in the quark dis- 
tributions. We cannot determine a priori the size of the polarized gluon distribution in a 



proton at a given Q^ value. The evolution equations for the polarized distributions, indicate 
that the polarized gluon distribution increases with Q^ and that its evolution is directly 
related to the behavior of the orbital angular momentum, since the polarized quark distri- 
butions do not evolve in Q^ in leading order. [^ Thus, one can assume a particular form 
for the polarized gluon distribution at a given Q\ and evolve it to the higher Q^ values of 
the data. Until we can experimentally check its consistency with data which are sensitive to 
AG'(a;, Q^) over a particular Q^ range, we must assume models for AG to analyze the spin 
properties of parton distributions. Initial analyses of the EMC data led to speculation that 
the integrated gluon distribution may be quite large, even at the relatively small value of 
g2 = 10.7 GeV2. 



Recent data from the E704 group at Fermilab [^ indicate that the polarized gluon 
distribution is likely not very large at the Q^ values of present data. With this in mind, 
there are two feasible models for a small AG, namely: 

(1) AG(x) =xG(a;), 

(2) AG(x) = 0. (62) 

The first implies that the spin carried by gluon is the same as its momentum, motivated 
by both simple PQCD constraints and the form of the splitting functions for the polarized 
evolution equations. The second provides an extreme value for determining limits on the 
values of the polarized sea distribution, assuming a positively polarized gluon distribution. 

There are models of AG which result in at least some negativity to the polarized 



gluons. Jaffe's model [^ is based on a constituent quark picture, where interactions with 



gluons are considered. This leads to a sizable negative polarization of the gluons and is 



naturally quite speculative. The model of Kochelev [^] is a non-perturbative model, which 
analyses vacuum fluctuations in the gluon field, called "instantons" . The kinematic analysis 
results in a polarized gluon distribution which has a negative component at small-x and a 
small positive component at larger x. The total integrated distribution is slightly negative, 
however. Both of these negative AG models would further compound the "spin crisis" prob- 
lem if the chiral invariant factorization scheme is used. 



The various gluon scenarios have been summarized by Di Salvo. |75| Nowak, et. al, 



have done an analysis on the data, including the axial anomaly, along with instanton-based 
q-q- interactions. |T^ They conclude that present data imply a positive AG and that the 
structure function g^ is highly sensitive to the sign of the polarized gluon distribution. This 
provides a possible test of these models. This is discussed in more detail in section III. 



2.3.2 Gluon Anomaly 

As we discussed in the factorization section, if a chiral-invariant factorization scheme is used, 
the polarized gluon distribution has an effect on the quark spin distributions via the axial 
anomaly, to be highlighted here. Using the same helicity notation as with the quarks and 
imposing parity invariance of the strong interactions we write 

AG(x, Q2) = G+/+(x, g2) _ G_/+(x, g2) (63) 

= G_/_(x,g2)_G+/_(x,g2), (64) 

which appears in the calculation of spin-related observables involving polarized protons in 
the same manner as the Agj(x, g^). The model of AG that is used has a direct effect on 
the measured value of the quark distributions through the gluon axial anomaly. |7^, ^ In 



QCD, the U(l) axial current matrix element A^ is not strictly conserved, even with massless 

quarks. Hence, at two loop order, the triangle diagram between two gluons generates a Q^ 
dependent gluonic contribution to the measured polarized quark distributions. This term 
has the general form: 

r(g2) = ^ll^^^ f AG{x, Q') dx, (65) 

ZTT JO 

where Nf is the number of quark flavors. Thus, for each flavor of quark appearing in the dis- 
tributions, the measured polarization distribution is modified by a factor: (Ag,;) —T{Q'^)/Nf. 
The quark spin contributions the depend indirectly on AG if the chiral invariant factoriza- 
tion scheme is used. In order for us to determine the quark contributions to the spin of the 
nucleons, it is necessary for us to know the relative size of the polarized gluon distribution. 
If we base our analysis solely on the naive quark model, then X) ^Q ~^ 1 a-iid AG may be 
quite large to be consistent with EMC data. If we consider the polarized distributions of Qiu 



et. al, [^ a reasonably sized AG is possible if the sea has a suitably negative polarization. 
In section III, we will consider two possible models for calculating the anomaly contribution: 



(1) AG* = xG (indicating that the spin carried by gluon is equal to its momentum) and (2) 
AG = 0, which is chosen to estimate a bound the distributions. Present data seem to imply 
that anomaly effects, and thus the overall integrated polarized gluon distribution, is limited 



at these energies. ||T0 



In the first moments of the parton densities, the contribution of polarized gluons to 
the integral of g^ is to produce an effective density 

(Ag,),., = (Ag,) - ^f^{AGiQ')) (66) 

for each flavor in the sea. The gluons change the measured net spin of the sea quarks and 
antiquarks by an amount 

E((Ag,)e.p - (Ag,)) = ^^^I^^{AG{Q')) (67) 

which may or may not be large. For a more detailed discussion of the anomaly term, see 



Cheng. [45 



2.3.3 Role of Orbital Angular Momentum in Nucleon Spin Content 

It is clear from the Jz = | sum rule that the angular momentum accounts for the amount 
of nucleon spin which is not carried by either the quarks or gluons. However, the orbital 
motion may play an even more important role than this implies. [^ The extended nature 



of the relativistic proton and its orbital motion may be responsible for the single spin asym- 



metries seen in the data ||T0[ for inclusive pion production. Furthermore, there are striking 
similarities between the inclusive hyperon polarization from unpolarized pp scattering and 
the orbital effects which explain the single spin asymmetries. This is a topic for further 



study and is being carried out by the Berlin group (Boros, et. al.). Troshin and Tyurin [79 
have mentioned that single spin asymmetries, which could be measured at the Relativistic 
Heavy Ion Collider (RHIC) at Brookhaven, would provide a good measure of the relative 
size of the orbital motion of the constituents. 



Ji, et. al., |80] have derived evolution equations for the quark and gluon orbital angu- 



lar momenta and have concluded that the asymptotic value for fractions of spin carried by 



quarks and gluons are 3ri//(16 + 3n/) and 16/(16 + 3?2/) as Q^ -^ oo. Ji |^ has also derived 
a sum rule for the orbital angular momenta and proposes a possible way of measuring these 
quantities. (See also, Radyushkin [0]). Thus, the orbital motion may be more interesting 
than originally thought. Naturally, if future experiments allow precise measurement of the 
spins carried by quarks and gluons, the Jz = ^ sum rule will provide another test of these 
angular momentum models. 



2.4 Models of the Polarized Sea 
2.4.1 Overview of Models 



There have been many approaches to extraction of the polarized sea from DIS data. |£2|, |8^, 
^ ^ Most of these are similar in the global analysis, where the goal is to find the fraction 
of nucleon spin carried by each flavor of the sea. Extraction of the x-dependent distributions 
generally follows along different lines, depending on the treatment of the polarized glue and 
the small-x behavior. 



The polarization of the sea occurs by gluons that are emitted by gluon Bremsstrahlung 
and by quark-antiquark pair creation. The total sea for three flavors is merely the sum of 
contributions from the flavors. This is written for the spin averaged and spin weighted cases 
as where 

S{x) = Us{x) + u{x) + ds{x) + d{x) + s{x) + s{x) (68) 

AS{x) = Aus{x) + Au{x) + Ads{x) + Ad{x) + As{x) + As(x), (69) 

where the valence quark contributions are omitted. The data give information about the 
integrated distributions which in the polarized case are 

{AS{Q^)) = {[Au,{Q^) + Au{Q^) + Ads{Q^) + Ad{Q^) + As{Q^) + As{Q^)]). (70) 

For the unpolarized structure functions, the strange quarks in the sea are often treated 
separately [0 to account for the excess of d over -u, which violates the Gottfried sum rule. 



In the polarized case, the heavier mass of the strange quarks will likely make them harder 
to polarize. In fact, the recent analyses of the DIS data agree that As is smaller than the 



lighter flavors, regardless of which approach is used to model the distributions. The details 
of extracting the integrated spin-weighted distributions from the data will be given in the 
next section. Here we will outline two possible ways of parametrizing the polarized sea quark 
distributions. 



In the approach of Goshtasbpour and Ramsey |4^, is is assumed that the power of 



(1 — x) is the same for both the polarized and unpolarized sea, indicating the same large x 
asymptotic behavior. For each flavor they assume 

Aqi{x) = r]ixqi{x), (71) 

where the polarization terms, rji, are a set of flavor dependent parameters, to be determined 
by data. Then, knowledge of the unpolarized distributions will yield the x-dependent po- 
larized densities. Then, the appropriate spin observables can be calculated to compare with 
data. In the Gehrmann/Stirling [^ and Bartelski/Tatur |2^ approaches, the polarized 
distributions are written as: 

S{x) = Ausx"^' (1 - xf^' (1 + ax^ + bx) (72) 

AS{x) = A,a;"==°(l - a;)^^^°(l + ax^ + bx), (73) 

and the parameters are extracted from data. 

In either of these scenarios, it is found that data imply a negatively polarized sea. 
This negative polarization can be generated in two possible ways: 

(1) the perturbative mechanism of the gluon axial anomaly, which breaks chiral symme- 
try, or 

(2) the non-perturbative instanton mechanism, which causes quark helicity flipping to induce 
the negative sea. 



Other theoretical arguments implying a negatively polarized sea involve non-perturbative 
spin-spin correlations and the large- A'^c chiral Lagrangian motivated by a Skyrme model of 
the proton. |7ll 



A detailed discussion of the sea distributions extracted from data is in the next major 



section. 



2.4.2 Small-,T Behavior 



The small-x behavior of the polarized distributions is crucial to both understanding the 
role of perturbative QCD in DIS and extracting the flavor dependent polarization densities 
from the data. Recent HERA data for the unpolarized distributions indicates growth of F2 
at small-x. [^] Since the polarized distributions are extracted from F2, this has a direct 



affect on the polarization of the proton in this kinematic region. (For a treatment of both 
unpolarized and polarized structure functions at small-x, see Webber [p| ) Historically, Regge 
theory predicted that f^f'" ~ x~", where the axial-vector meson trajectory, a was in the range 
—0.5 < a < 0.0. Data indicate that this prediction is likely valid for Q^ < 1 GeV^, but that 
perturbative effects are more appropriate for higher Q^. 

There are two popular approaches to analyzing the small-x behavior of the polar- 
ized structure functions. The first is a standard resummation of the GLAP Q^ logarithms: 
En,n^a.(i)"(lnQ')"^, implying H 

(in I i 1)" < (^r < ^''. (74) 

for some positive p and q. The second is BFKL inspired resummation of (1/x) logarithms: 



m 



I]n,m'^s(^o)"(ln(l/x))™ at fixed Q^ where the non-singlet and singlet structure functions 
behave as: IBSl 



ar~'=-°\'^\r (75) 

The resummation of higher order corrections is important in understanding the small-x be- 
havior of the polarized structure functions. ||6^ Preliminary data from the SLAG E154 
experiment can be fit to a power law g^ ~ x~°'^ for 0.02 < x < 0.1, but Ratcliffe has pointed 
out that the data can be fit equally well to a form: gf" ~ x~°'^(l —Ax). ||8^ This is consistent 



with an isospin decomposition of data done by Soffer and Teraev, [R^ who give a small-x 



power of x~^''^^. The small-x extrapolation can give a net value of the integral /g dxQi which 
differs by up to a factor of two. This makes a large difference in the extraction of spin infor- 
mation from this data. So far, HERA data has not been able to differentiate between these 
scaling models for F2, but future experiments at HERA or the LHC could reach lower values 



of X at high enough Q . |88| A recent analysis of gi via an all-order resumming of the 



0(0;^"''^ In x) terms in the singlet evolution, indicates a large uncertainty in the behavior in 
these structure functions at small-x due to uncalculated correction terms. Thus, it appears 
that both theoretical and experimental work must be done in order to isolate the small-x 
behavior of the polarized structure functions, and hence, the spin distributions. 
A crucial problem here is extrapolation of the structure function gi{x) to x — i> 0. The region 
of small X is particular interesting since it provides an insight into interface of the perturba- 
tive and nonperturbative regions of QCD. The small-x behavior of structure function gi{x) 
has been described traditionally in the Regge model by the contribution of ai trajectory with 
intercept ^^^(O) ~ 0. However, the SMC data point out that gi{x) might increase at small x. 
Theoretical background for the Regge extrapolation is also questionable, since we deal with 
the amplitudes with virtual external particles. Indeed perturbative QCD evolution gives 



another form for gi at small x, i.e. gi{x) ~ exp win 1/x. Other forms of this dependence are 
allowed in general Regge analysis with account for cut-contributions, two-gluon model for 
Pomeron. Even strongly rising at x — ;> dependencies such as gi ~ In^ x/x are possible. In 
the latter case the integral for the first moment of gi is divergent. An important question 
here is the role of unitarity for the amplitudes with virtual external particles, i.e. whether it 
provides any restrictions to the growth of gi. Thus the problem of extrapolation to x — i> is 
important. The only way to resolve this problem is with the experimental measurements in 
the region of x ~ 10~^ — 10~^. Such measurements are possible with the polarized proton 
beam at HERA. 



Another potential contribution to the spin analysis at small-x is the possible Q^ scal- 
ing violation of the asymmetries A[ (i=p,n,D). A detailed analysis done by Cluck, et. ai, 
B9| indicates that this slight scaling violation could exist at most x < 0.25 for Q^ < 4 CeV^. 



The potential problem here is that the experimental analysis is done using the assumption 
that the asymmetries are independent of Q^. It is unclear that this effect would be signifi- 



cant, since the scaling violation is small and the values of g{ are more sensitive to the small-x 
fit than the scaling of Ai. This will be discussed in section III. 

All sets of data are limited in the range of Bjorken x and thus, the integrals must 
be extrapolated to a; — > 0. Thus, the possibility of existance of a Regge type singularity 
at X —> is not guaranteed in the analyses. A significant singularity could raise the value 
of Qi towards the naive quark model value and could account for some of the discrepancy 
between the original EMC data and the Ellis- Jaffe sum rule. In light of the recent HERA 
data, there is the possibility that the increase in F2 at small x, even at the lower Q^ values 
of the E142/E143 data, could indicate a change in the extrapolated values of these integrals. 
These possibilities are a topic for future study. It is also possible that the overall effect of 
F2 on Qi will not alter the integral by any more than the present experimental errors. The 
shape of the polarized gluon distribution at small-x affects the anomaly term, and thus the 
overall quark contributions to the integrals. Future experiments can shed light on the size 
of this effect, a detail discussed later. There is still controversy as to whether present data 
show that anomaly effects are limited or not, thus leaving open the question of the size of 
the polarized gluon distribution at these energies. This will be discussed in detail in the next 
section. 



3 EXPERIMENTAL CONSEQUENCES 

3.1 Experimental Overview 

3.1.1 Recent Data from SLAC, CERN and DESY 

The most recent generation of Deep-Inelastic Scattering (DIS) experiments began with the 
European Muon Collaboration (EMC) at CERN in 1987. [Q Their measurement oi Qi yielded 
results which were in disagreement with the Ellis- Jaffe sum rule (EJSR). Naturally, this mo- 
tivated many theoretical analyses, as well as plans for further polarized DIS experiments to 
test the models which were devised to explain the discrepancy. As a result, experimental 
groups at CERN and SLAC performed a number of these experiments. Since 1992, these 
groups have succeeded in measuring the corresponding neutron and deuteron structure func- 
tions as well as measuring Qi more precisely than the EMC. In all cases, the statistical 



and systematic errors were significantly decreased and a wider range of Bjorken-x values 
was probed. This improvement has been promoted by both technological developments and 
increased running time. 

Beginning in 1993, the SMC group measured gf, while the E142 experiment measured 
(?". This provided a complementary set of measurements to test both the Bjorken sum rule 
(BSR), which measures the difference of their integrals, and the Ellis- Jaffe sum rule, (EJSR), 
which predicted their separate integrals. In 1994 and 1995, both the SMC and E143 groups 
measured Qi and gf to check consistency with former results and provide a "world" average 
for these quantities, so that a more thorough theoretical analysis could be carried out. 
A number of theoretical teams performed these analyses. [^, ^, ^ ^ ^ Before the 
recent E154 and HERMES data were released, the neutron data seemed to yield different 
implications for the polarized strange sea than the proton and deuteron data. These groups 
have measured g^ to higher precision and over a wider range of x so that a revised value was 
given. This has brought the neutron data closer to the implications of the other data, but 
this issue is still not resolved. 

A comparison of experimental results and their corresponding theoretical analyses 
will be summarized in this section. The physical implications of the analyses will then be 
discussed, along with a planned set of future experiments which will help to resolve some of 
the unanswered physics questions. 

A summary of key measurements in polarized DIS is given in Table I. Note the com- 
plementary measurement of the structure functions as well as the increased coverage of the 



X range and decreasing experimental errors. Note that the E154 ||63| and HERMES [§[ data 
are preliminary. 



Table I: Recent DIS Experimental Data 


Experiment 


Target 


iQ') 


x-range 


Ttarget 


EMC (87) 


P 


10.7 


0.1 <x< 0.7 


0.126 ±0.010 ±0.015 


SMC (93) 


d 


4.6 


0.006 < X < 0.7 


0.023 ±0.020 ±0.015 


SMC (94) 


P 


10.0 


0.003 < X < 0.7 


0.136 ±0.011 ±0.011 


SMC (95) 


d 


10.0 


0.003 < X < 0.7 


0.034 ±0.009 ±0.006 


E142 (93) 


n 


2.0 


0.03 < X < 0.6 


-0.022 ±0.007 ±0.006 


E143 (95) 


P 


3.0 


0.03 < X < 0.8 


0.127 ±0.004 ±0.010 


E143 (95) 


d 


3.0 


0.03 < X < 0.8 


0.042 ± 0.003 ± 0.004 


E154 (96) 


n 


4.5 


0.014 < X < 0.9 


-0.037 ±0.004 ±0.010 


HERMES (96) 


n 


3.0 


— 


-0.032 ±0.013 ±0.017 



Most of the CERN experiments have used polarized muons, while the SLAC and 
DESY (HERMES) experiments have used electrons, being primarily electron accelerators. 
As Frois pointed out, |^| electron data tend to be more accurate due to the higher intensities 
possible, but provide a less broad kinematic range than muons, due to their smaller energies. 
The systematic errors in each are similar. Thus, these experiments provide complementary 
information about the polarized structure functions. 



The targets typically used for these experiments are: ammonia crystals (NH3) for the 
proton data, ^He for the neutron data and deuterated butanol (^^NDs) for the deuteron data. 
More detailed information on targets is contained in workshop proceedings: D.G. Crabb, N. 
Horikawa and V.G. Luppov in the VI Workshop on High Energy Spin Physics, Protvino 
(Vol. 2), 1996, and S. Goertz, V.G. Luppov, B. Owen in the 12th International Symposium 
on High Energy Spin Physics, Amsterdam, September, 1996. See also the reference of Crabb 
and Day. |9^ Szwed |9^ has done a partial study on the effect of the nuclear targets on the 
reliability and understanding of the data. 

Other recent technical developments in polarization experiments include improve- 
ments on Siberian snakes used to retain polarization of accelerating protons ||9^, better 
polarized ion sources and detector improvements. pSl 



3.2 Extracting Results from Data 

3.2.1 Extraction of Polarized Sea and Gluons from Data and Sum Rules 



Being the earliest analysis, the Close and Roberts approach |^ used leading order structure 



functions and used the earliest data (EMC proton, E142 neutron and SMC deuteron). They 
start with the BSR and use average values for the F and D weak decay constants, discussed 
in section 11. They introduce an abitrary parameter for the possible higher-twist corrections. 
A value for the total spin carried by quarks, Aqtotai, was found for each experiment. The 
resulting error bars were quite large, since data were over a smaller kinematic range and the 
experimental errors were larger. Further, the higher order QCD and higher-twist corrections 
were not known as well as at present. 



Ellis and Karliner [^ expanded this analysis by including higher order QCD and 
higher-twist corrections and taking a world average of the data available as of their writing. 
This included the EMC, E142 and E143 experiments as well as the SMC proton and deuteron 
experiments. They used the sum rules to extract flavor dependent information about the 
proton spin. 

Cheng and Li [^] took a more theoretically motivated approach by considering a 
chiral quark model with the Gottfried sum rule violation built in (asymmetry of anti-down 
to anti-up quarks). Their experimental inputs were the value of this asymmetry, measured 
by the NMC group at CERN ||3^ and the (Ttttv factor from pion-nucleon experiments. Their 
polarized sea was completely SU(3) symmetric. 



The approach of Goshtasbpour and Ramsey PS[ included higher order corrections to 



the structure functions and did a complete flavor dependent analysis, including two different 
gluon models to investigate the effects of the gluon anomaly. They also break the SU(3) 
symmetry of the flavor dependent sea to explicitly separate the strange sea contribution. 
This was repeated for each of the experiments separately and therefore indicated where the 
data were consistent and where they disagreed in terms of physical implications. Being a 
more explicit analysis, it reflects many of the techniques used by the theoretical groups to 
analyse the data. It will therefore be outlined shortly as an example. Their 1996 update 
includes preliminary data from the E154 and HERMES experimental groups, which will be 
covered later. 



The Ball, Forte and Ridolfi analysis [^] places emphasis on the low-x behavior of the 
structure functions and does a moment analysis of their evolution. They assume complete 
SU(3) symmetry and extract both quark and gluon spin information from the data. A 
comparison of the results from these theoretical approaches will be discussed later. 

To give an indication of how the data and sum rules are used to extract the spin 
information about the quark and constituents, we will explain the basic approach of Gosh- 
tasbpour and Ramsey (G-R). The basic source of information for fixing the parameters in 
the models is the data on the longitudinal spin-spin asymmetry in deep-inelastic lepton- 
proton scattering and the sum rules discussed in section II. The measured integral of gi in 
each experiment constrains the appropriate spin parameters. The higher-twist corrections 
53| appear to be negligible at the Q^ values of the data, so they are not included. The 



additional constraints are provided by the axial-vector current operators, 03, og and oq, as 
discussed earlier. The BSR is used to extract an effective P = J^ gi{x) dx from all data. 

The model of AG that is used has an affect on the quark distributions through the 
gluon axial anomaly, which was discussed earlier. In Table II, two models for AG are con- 
sidered: (1) AG = xG and (2) AG = 0. The E154 and HERMES data are preliminary, as 



reported in the Amsterdam symposium. p3 



The key elements of the G-R approach are: 

• determine the valence contribution to the spin using the BSR 

• find sea integrated parton distributions for each flavor by breaking the SU(6) symmetry 

with the strange quarks and using the sum rules with data as input 

• include higher order QCD corrections and the gluon anomaly for each flavor 

• discuss similarities and differences between the phenomenological implications of the 
different experimental results, and 

• suggest a set of experiments which would distinguish the quark and gluon contributions 
to the proton spin. 

This approach differs from that of others in that the sum rules are used in conjunction 
with a single experimental result to extract the spin information and the flavor symmetric sea 
is broken while anomaly contributions are included via the gluon models. The sea breaking 
parameter, e, is defined by: 

/^Usea = Au = Adsea = Ad = {1 + e)As = (1 + e) As. (76) 



The analysis (for each polarized gluon model) proceeds as follows: 

• Extract a value of P from either the data directly or via the BSR in the form of 
equation (77), 

• use Eqn. (44) to extract a^. Then the overall contribution to the quark spin is found 

from (Aqtot) = ^o + T, where F is the gluon anomaly term in equation (65), 

• use the value Og from the hyperon data to extract As for the strange sea, 

• find the total contribution from the sea from {Aqtot) = (Aqy) + (AS), 

• determine the SU(3) breaking parameter, e and the distributions {Au)sea = {^d)sea 
from equation (76) and the strange sea results and 

• finally, extract L^ from the Jz= 1/2 sum rule. 

3.2.2 Comparison of Results from Different Experiments 



Data from SMC [g, SLAC [17, 63] and DESY [H are used to extract information about the 
flavor dependence of the sea contributions to nucleon spin. We can write the integrals of the 
polarized structure functions, P = /g g\ dx in the terms of the axial- vector currents as: 



P 



jn 



jd 



g{{x)dx 



g1{x)dx 



;i - I^d) 



12 36 9 



(1 



a. 



(77) 



' 12 36 9 



a„ 



gf{x)dx 



36 9 



(1 



a„ 



;i - Iu^d) 



where ujd is the probability that the deuteron will be in a D-state. Using N-N potential 
calculations, the value of Ur) is about 0.058. [Q The BSR can then be used to extract an 
effective P value from all data. Comparison of the /^ r r values from each experiment gives a 
measure of the validity of the BSR. 

Since the evolution splitting functions for the polarized distributions have an addi- 
tional factor of X compared to the unpolarized case, early treatments of the spin distributions 
assumed a form of: Aq{x) = xq{x) for all flavors. This form of the distributions has been 



(Aqsea)e 



^ for 



compared to those extracted from the recent data, using the defined ratio ri = , \ 

each flavor. Any deviation from rj = 1 would indicate that the early models for generating 
the polarized distributions are inaccurate. The results are given in Table II. 



Table II: Integrated Polarized Distributions: 
AG = xG (above line), AG = (below line) 



Quantity SMC{P) SMC{r 



< Au >sea 

< As> 

< Au >tot 

< Ad >tot 

< As >tot 
Vu = Vd 

Vs 

e 

r 

p 

< Aq >tot 
<AG> 



-.077 
-.037 
0.85 
-.42 
-.07 
-2.4 
-2.0 
1.09 
0.06 
0.136 
0.36 
0.46 
-.14 



-.089 
-.048 
0.82 
-.43 
-.10 
-2.8 
-3.0 
0.84 
0.06 
0.129 
0.29 
0.46 
-.11 



E154(J" 

-.063 
-.020 
0.87 
-.39 
-.04 
-1.9 
-1.2 
2.10 
0.08 
0.134 
0.45 
0.45 
-.18 



^143(1'=' 

-.068 
-.028 
0.87 
-.40 
-.06 
-2.1 
-1.6 
1.41 
0.08 
0.131 
0.41 
0.44 
-.15 



HERMES 

in 

-.050 
-.010 

0.90 

-.36 

-.02 

-1.5 

-0.6 

4.00 

0.07 
0.135 

0.52 

0.44 

-.22 



< Au >tot 

< Ad >tot 

< As >tot 

P 

< Aq >tot 

r 



.83 
-.44 
-.09 
.136 
0.30 
0.00 
0.35 



.80 
-.45 
-.12 
.129 
0.23 
0.00 
0.39 



.85 
-.41 
-.07 
.134 
0.37 
0.00 
0.32 



.84 
-.43 
-.08 
.131 
0.33 
0.00 
0.35 



-.39 
-.04 
.135 
0.45 
0.00 
0.28 



3.3 Consequences of the Results 



3.3.1 Physics Consequences 

From these resuhs, it is obvious that the naive quark model is not sufficient to explain the 
proton's spin characteristics. Nor is the simple model for extracting the polarized distribu- 
tions accurate. Some conclusions which can be drawn from the data are: 

(1) The total quark contribution to proton spin is between 1/4 and 1/2. The errors in 
generating these results are due mostly to experimental errors and determination of which 
model of the polarized gluons to use. 

(2) The up and down sea contributions seem to agree within a few percent. However, 

the proton and deuteron data imply a larger polarized sea with the strange sea polarized 

greater than the positivity bound. ||58| Interestingly, the SMC proton data are consistent 

with a recent lattice QCD calculation of these parameters. |7^ The results from these data 

can be categorized into distinct models, characterized by the size of the non-zero polarized 
sea. 

(3) The values of ?7 deviate considerably from unity for most of the data, implying that 
the relation between unpolarized and polarized distributions is more complex than originally 
thought. 

(4) This analysis implies that the anomaly correction is not large. If the anomaly 
term were larger, due to a large AG, the strange sea would be positively polarized, while the 
other flavors are negatively polarized. There is no known mechanism that would allow this 
cross polarization of different flavors. These data imply that AG is of small to moderate size. 
Further, even if there are higher twist corrections to the anomaly at small Q^, the anomaly 
will not reconcile differences in the flavor dependence of the polarized sea. 

(5) The orbital angular momentum extracted from data is also much smaller than 
earlier values obtained from EMC data. In fact, a small AG model implies a correspond- 
ingly small orbital angular momentum, although its sign is still in question. 

(6) The extracted P value is comparable for all data and well within the experimental 
uncertainties. This implies agreement about the validity of the Bjorken Sum Rule. This has 



been done here by using the BSR to extract an effective P, in contrast to other analyses, 
which use data to extract the BSR. There is general agreement that the BSR (and thus 
QCD) is in tact. 

Clearly, these experiments have contributed to the progress of understanding the rel- 
ative contributions of the constituents to the proton spin. They have probed to smaller x 
values, while decreasing the statistical and systematic errors. This, coupled with theoretical 
progress in calculating higher order QCD and higher twist corrections have allowed us to 
narrow the range of these spin contributions. Although the flavor contributions to the proton 
spin cannot be extracted precisely, the range of possibilities has been substantially decreased 
(see Table III). The main differences are the questions of the strange sea spin content and 
the size of the polarized gluon distribution. Obviously, more experiments must be performed 
to determine the relative contributions from gluons and various flavors of the sea. 

Table III: Ranges of Constituent Contributions to Proton Spin 



Quantity 


EMC results 


Post-SMC/SLAC 


(A«).ea 


-0.077 


-0.089 


(As) 


-0.037 


-0.028 


{^^)tot 


0.85 


0.80 ^ 0.90 


{Mtot 


-0.42 


-0.35 ^ -0.45 


(A3),,, 


-0.25 -^ 


-0.12^0 


P 


0.126 


0.136 


(Mtot 


0^ 1 


0.2 ^ 0.5 


(AG) 


0^6 


-^ 1.50 


L, 


0^6 


-^ 1.25 



Table IV shows a comparison of various recent approaches to extraction of the spin in- 
formation from data. The key is as follows: BFR [||]; CL ||; CR pl; EK || and GR. H 



Table IV. Comparison of Results with Different Models 



Quantity /Model— 


BFR{9Q) 


CL{95) 


CR{93) 


EK{95) 


GR{9Q) 


< Am >tot 


0.88 


0.79 





0.83 


0.86 


< Arf >tot 


-0.38 


-0.32 





-0.43 


-0.40 


< As >tot 


0.00 


-0.10 


— 


-0.10 


-0.06 


P 


0.122 





0.126 


0.133 


0.133 


< Ag >tot 


0.50 


0.32 


0.38 


0.30 


0.40 


<AG> 


1.50 


0.00 








0.45 


L, 


-1.25 


0.34 








-0.15 



It is clear from these results, that even with varied approaches and assumptions, the 
up and down polarized distributions are fairly consistent. A slightly different approach was 
taken by Gliick, et. ai, |^ in a next-to- leading order analysis. Their results were com- 
parable to BFR, except for the polarized strange sea, which agreed with the GR approach. 
The range of As and AG are quite considerable in these models. Part of the problem is the 
small-x contributions to the structure functions, which imply different gluon contributions. 
Also the lack of knowledge of the orbital component of motion prevents us from making a 
statement about the glue or the strange sea. 



3.3.2 a;-Dependent Distributions 



In order to generate the x-dependent distributions, there are two approaches which were 
mentioned in section II. Three groups have extracted these distributions directly from the 



data. ||97| , |23| , |27[| All give good agreement with data, but they differ significantly in the 
small-x region. Some are consistent with Regge behavior, while others are not. Goshtas- 
bpour and Ramsey |^, 48 1 have used the unpolarized distributions with their extracted value 
of 7] and the assumption that: Aq{x) = rixq{x) for each of the sea flavors. For the valence 



distributions, they have used the model of Qiu, et. ai. |2^ There is no reason a priori 
to suspect that a global fit to the integrated distributions should imply a satisfactory x- 
dependent fit to the data. However, figs. 1 through 3 indicate that this form gives very 
good x-dependent parametrizations for the polarized distributions, consistent with data and 
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Figure 2: The x dependent proton and neutron structure functions, g\ and g^ generated 
from Eqn. (71) and data. The dotted hue represents the GRV-generated distributions and 
the sohd hne, MRS-generated distributions. 

Regge behavior. 



The differences between the these sets of distributions are at small-x, where the data 
is most uncertain. It is clear that more DIS experiments should be performed to probe 
very small-x to distinguish between models and to address the controversy regarding which 
contributions to g\ dominate in this kinematic region. It will be required, however, that the 
error have to be minimized to distiguish between the various possible powers of small-x for 
the polarized structure functions. 



3.3.3 Open Physics Questions 



The proceeding analyses emphasize a number of key physics questions which the latest DIS 
experiments have raised. These include, but are not limited to: 

• the size and sign of the polarized gluon distribution 

• the amout of spin carried by the strange sea 
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Figure 3: The x dependent deuteron structure function, gf, using the same technique as in 
figure 2. 



the light quark flavor dependence of polarization; is there an analogy to violation of 

Gottfried sum rule in Am and Arf? 

the connection, if any, between the unpolarized and the polarized distributions 

the role of the higher order corrections at low x and Q^-both perturbative and higher 

twist (non-perturbative) 

the role of the orbital motion, L^; does it agree with the expected asymptotic values? 
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Theoretical models which agree with present data still disagree on many of these 
points. It is therefore up to the experiments to decide which of these explain the various 
aspects of constituent spin contributions to the nucleons. In the next section, we outline the 
various polarization experiments which will address these issues. 



3.4 Future Experiments 
3.4.1 Introduction 

There are a number of experiments which are technologically feasible that would supply some 
of the missing information about these distributions. Detailed summaries can be found in 
references by Ramsey |^ and Nurushev. ||^ We have seen that the existing data have 



enabled us to formulate appropriate questions which probe the spin properties of nucleons. 
However, there are a number of questions which remain unanswered and will only be acces- 
sible with more data at different energies and momentum transfers. Fortunately, there are 
a number of experimental groups that are planning polarized beam experiments at exist- 
ing accelerators. With recent advances in polarized beam, target, detector and accelerator 
technology, it is now possible to do these experiments at higher energies and momenta in 
order to study the physics over a large kinematic range. The large average luminosities of 
these experiments and the success of Siberian Snakes makes all of the following feasible. 
This section will include some of the proposed experiments related to the spin structure of 
nucleons, in light of the questions presented in the last section. 

One of the ways to categorize polarization data is as follows: 
• deep-inelastic scattering of polarized leptons (e, /x) on polarized nucleon targets (p, n, 
d) 



• photo-production of jets in high energy polarized ep coUiders 

• production of pions and direct photons in polarized pp and pp scattering 

• charmed meson production {J/ip and x) in PP collisions 

• direct photon and jet production in polarized pp collisions 

• lepton pair production (Drell-Yan) in polarized processes 

• heavy baryon (hyperon) production in unpolarized pp collision 

All of these are designed with the measurement of particular distributions in mind. 
Some can also provide crucial tests of QCD. There are a number of these presently planned at 
the following locations (alphabetically): (1) CERN (Switzerland), (2) DESY (Germany), (3) 
LISS (Indiana, USA), (4) RHIC (Brookhaven, USA), (5) Serpukhov (Russia) and (6) SLAG 
(USA). A partial list of experiments with their corresponding energies is given in Table V. 
This is to give a general idea of the wide range of energies and kinematic regions to be covered. 



Table V 



Experiment 


Location 


Energy ^/^{GeV/c) 


Luminosity (cm ^ s ^) 


HERMES 


HERA (DESY) 


30 (e) on 820 (p) 


2 ■ 1031 


SPIN 


HERA (DESY) 


820 


2 ■ 10^1 


RHIG 


Brookhaven 


60 -^ 500 


2 ■ 10^2 


LISS 


Indiana 


20 


1 ■ 10^2 


GOMPASS 


GERN 


120 


1 ■ 10^2 


NEPTUN-A 


Serpukhov 


400 


1 ■ 10^1 


E155 


SLAG 


48(e) 


— 



In the following discussion, we will describe some of these experiments with regard to 
the physics that they probe, namely: 

• a. Ag measurements (valence and sea) 

• b. AG measurement 

• c. Lz determination 

• d. Nuclear measurements - F^, g^, hi ' 

• e. higher order corrections and other tests of QGD 



A partial list of the experiments for these categories is in Table VI. These will be 
discussed in the following subsection. The experiments discussed here represent a sampling 
of those which directly relate to the subject of the spin structure of nucleons. Many of these 
are just a small fraction of the polarization experiments which can be performed at these 
accelerators, but they form an integral part of the program. 

Table VI 



Experiment 


Proposed Type 


Measured Quantities 


Distribution 


HERMES 


DIS 


^l9{ 


Ag, AG 


E155 


DIS 


glgf 


Ag, A(? 


SPIN 


Elastic pp 


An-, Ann 


Helicity NG 


RHIC/COMPASS 


Gharm prod. 


A" 


AG 


RHIC (STAR) 


Jet, vr, 7 prod 


A(Tl, All 


AG 


RHIC (PHENIX) 


Drell-Yan 


^DY 


A^ 


LISS 


Inelastic 


0"L, 0"T, ^<^L 


AG 


SLAG 


Gharm prod. 


A" 


AG 


SLAG 


W± prod 


j^W 


Agi 


LHG 


W± prod 


J^W 


AG 



3.4.2 Valence and Sea sensitive experiments 

Deep Inelastic Scattering: The E155 experiment has been approved at SLAG. This experi- 
ment is designed to probe slightly smaller x while greatly improving statistics and system- 
atical errors. With lower error bars at small x, the extrapolation should achieve a more 
accurate value for the integrated distributions and narrow the ranges of constituent spin 
contributions even further. 



There has been considerable discussion about performing the GOMPASS polarization 
experiments at the LHG at GERN. Depending on the approved experiments, there is the 
possibility of probing small x and doing polarized inclusive experiments to measure both sea 
and gluon contributions to proton spin. These could be made in complementary kinematic 
regions to those of the other accelerators. There are tentative plans to do polarized W^ 



production, which provides a measure of the x-dependent sea distributions. Polarized W^ 
production is also planned at SLAC and would provide useful flavor dependent sea informa- 
tion in a slightly different kinematic region than that of CERN. 

The latest experiments at HERA in Hamburg have accelerated a large flux of polar- 
ized electrons from the storage ring and collided them with a gaseous target. The gaseous 
target has helped to eliminate some of the systematic errors characteristic of solid targets, 
which were used in the other experiments. This can be repeated for proton targets to test the 
BSR. With more events and the lower error bars at small x, it will be easier to extrapolate 
Qiix) and achieve a more accurate integrated value, {g{)- Thus, comparison to the Bjorken 
and Ellis- Jaffe sum rules will be more accurate, as will the ability to determine the polarized 
sea values from this data. p3|, |100| 



Anselmino, et. al, ||101|| have proposed doing charged current interactions [l^p -^ 
z/X), which could measure various linear combinations of flavor dependent polarized dis- 
tributions. These experiments are feasible and could put further constraints on the quark 
contributions to spin. 

Tests of the valence quark polarized distributions can be made, provided a suitable 
polarized antiproton beam of sufficient intensity could be developed. [[102|] This would pro- 
vide a good test of the Bjorken sum rule via measurement of (Ag^,) and the assumption of 
a flavor symmetric up and down sea. 

Lepton pair production (Drell-Yan) processes provide another clean measure of the 
polarized sea. ||103|| The Relativistic Heavy Ion Collider (RHIC) at Brookhaven is designed 



to be an accelerator of both light and heavy ions. ||104|| The high energy community has 
proposed that polarized pp and pp experiments be performed, due to the large energy and 
momentum transfer ranges which should be available. The energy range will be made in 
discrete steps between 50 and 500 GeV, and the momentum transfer range also covers a 
wide kinematic region. There are two main proposed detectors, STAR and PHENIX, which 
have different but complementing capabilities. Polarized Drell-Yan experiments are planned, 
which would give reasonable estimates to the polarized sea for each flavor. The PHENIX 
detector is suitable for lepton detection and the wide range of energies and momentum trans- 



fers could yield a wealth of Drell-Yan data over a wide kinematic range. The x-dependence 

of the polarized sea distributions could then be extracted to a fair degree of accuracy. Kamal 

105 has calculated next-to-leading order (NLO) corrections to Drell-Yan processes and has 



made predictions for RHIC energies. 

There has been a major effort to propose a double polarized pp mode at HERA, which 
would accelerate polarized protons. There is a wealth of both inclusive and exclusive ex- 
periments that can be done there at an energy range complementary to that of the other 
accelerators. This would require a major upgrade and installation of Siberian snakes, but 

the physics output potential is great. This initiative would be labeled HERA-A^. For a 

comprehensive review of the prospects at HERA-A^, see the proceedings of the workshop 
at Zeuthen ||106|| and references by Anselmino, et. al, [|107|| and Nowak. ||108|| Gehrmann 



and Stirling |2^ have calculated the NLO Drell-Yan asymmetries for both unpolarized and 

polarized experiments at HERA-A^ energies. Combining these with the RHIC measurements 
would give a good measure of the polarized sea over a large energy range. 



3.4.3 Gluon sensitive experiments 

Knowledge of the polarized gluon distribution is important for both the anomaly contribu- 
tion to nucleon spin and to its over-all spin contribution via the Jz sum rule discussed earlier. 
Thus, experiments sensitive to AG are a high proirity item. The following is a brief discus- 
sion of the possible gluon sensitive experiments and their corresponding proposed lacations. 

(1) Jet production in polarized ep collisions: It has been suggested that AG could be 
measured in ep collisions which produce one or two jets of hadrons from the photon-gluon 
fusion process. ||109| , |110|| This experiment could be done at HERA, but the experted asym- 



metries, even for a large polarized gluon distribution are only at the few percent level and 
may be difficult to measure with a sufficient degree of certainty to distiguish between the 
gluon models. 

(2) Jet production in pp collisions: Even before the spin crisis was popular, it was 



known that jet production in polarized pp collisions could be a sensitive measure of AG. [ |111 



The STAR detector at RHIC is suitable for inclusive reactions involving jet measurements. 
This would provide an excellent measurement of the Q^ dependence of AG due to the large 
range of energies available there. The NLO corrections to photon plus jet production have 
been calculated by Gordon. [p.l2|| The corresponding asymmetries range from a few percent 
for small gluon polarization to about 30% at HERA energies, if AG is large. Should DESY 
proceed with plans to polarize their proton beam, this experiment could be performed there, 
complementing the kinematic regions covered by RHIC and CERN. 

(3) Direct photon and double photon production: A clean signature for AG, but one 



that is harder to measure, is that of direct photon production, [pp ^ 7 + X) . p9| , 113 The 



NLO corrections have been calculated for HERA-X |n|, |TT5[ and RHIC. |TT6[ Although the 



asymmetries are only very large for the largest gluon models, the signal is clean and could 
distinguish if AG is large. Numerical simulations have been done in NLO for the double 
photon production processes {pp — > 77 + X) . [|117|| At RHIC energies, the asymmetries can 



be from a few percent to 20% for the larger AG, but the cross sections arequite small for the 
optimal momenta (picobarns or smaller). Thus, it is unclear whether this is viable unless 
very high luminosities can be reached. 

(4) Charm production in polarized collisions are also sensitive to AG and should be 

performed at both RHIC and HERA-X. The double spin assymetries for J/^ production 
have been calculated and are sensitive to the color-octet contribution. [ |118| , |1 19| | These asym- 
metries are only on the few to 10% level. The two-spin asymmetries for x production have 
also been calculated. ||120|| These asymmetries are quite small unless the polarized gluon 
distribution is quite large ((AG) ~ 6). Thus, x production could distinguish an extremely 
large distribution from other models. Open charm production from lepton-hadron scattering 
has been proposed at LHC by the COMPASS group. [^ 121 , 122 | This provides another 



method to measure the photon-gluon fusion process and has very small uncertainties asso- 
ciated with the proposed experiment. These asymmetries can be quite large and provide a 
good potential for extending our knowledge of AG. 

(5) Inclusive reactions involving pion production would be alternate tests of the Q^ 
dependence of AG. [|102|| These could be done at any of the aforementioned accelerators and 



would provide a good cross check of measurements of AG in similar kinematic regions. 



3.4.4 Experiments Probing Higher order QCD Effects 

There are other polarized experiments which would provide tests of QCD and give a measure 
of some of the higher twist effects which were previously discussed. 

(1) Elastic scattering, especially at a large —t range, could shed light on both helicity 
non-conservation at the hadronic level and non-perturbative long range effects. Elastic ex- 



periments have been proposed for HERA, RAMPEX at Serpukhov and LISS in Indiana. [ |123|| 



(2) Recently, a proposal for a new light ion accelerator was announced, which will 
specialize in polarization experiments. [[L24|| The Light Ion Spin Synchrotron (LISS) would 
be located in Indiana to perform a variety of polarization experiments for both high energy 
and nuclear physics. The energy range would be lower that most other experiments, thus 
complementing the kinematic areas covered. Furthermore, both proton and deuteron beams 
could be available to perform inclusive scattering experiments. They propose to measure 
longitudinal and transverse cross sections and spin asymmetries, which will address the nor- 
malization of the proton wave function. Elastic scattering measurements of Aj^ at moderate 
momentum transfer —t could give valuable information regarding helicity non-conservation 
in this region. 

(3) Measurement of the transverse spin and transversity distributions would probe 



higher- twist non-perturbative effects as well. |125| The higher twist parton distributions 



could also be measured at RHIC and HERA, pe], ^, [128 



(4) Single spin asymmetries in DIS (where only one of the scatterers is polarized) 
have been calculated by many groups. These provide helicity conservation tests as well as 



measures of higher twist processes. ||129| , |79| , [78[] These experiments could be done at HERA. 
The pp single spin asymmetries can also be done at HERA, RHIC and the LHC. These are 
good measurments of higher twist contributions to the structure functions. The pp single 
spin experiments are also planned by the RAMPEX collaboration (Russian-AMerican Po- 
larization Experiment). ||130|| 



This is a summary of some of the key experiments which will answer the questions 
posed throughout this review. There are many other important polarization experiments 
which test QCD and non-standard physics, but they will not be covered here. However, 
from this list alone, it is easily seen that polarization experiments are crucial to our under- 
standing of the most fundamental properties of the elementary particles. 



3.5 Conclusion 

This review has outlined the key theoretical and experimental elements which have con- 
tributed to our understanding of the constituents' contribution to the spin of nucleons. 
Analyses of existing data indicate that the spin structure of nucleons is non-trivial and has 
led to the formulation of a crucial set of questions to be answered about this structure. 
The key remaining questions are related to the strange sea and gluon polarizations. The 
experiments discussed here can be performed in order to shed light on these questions. In 
performing these experiments, an added benefit is that crucial tests of QCD and the quark 
model will be simultaneously be made. The questions of factorization and validity of the sum 
rules are among the questions which can be addressed. In the past few years, theorists have 
made considerable progress in calculating higher order corrections to the appropriate sum 
rules and structure functions and experimentalists have added more accurate and compre- 
hensive DIS data. This has contributed to narrowing the range of possible spin contributions 
from the constituents, but has also raised other questions about spin structure of nucleons. 

There is still much work to be done in constructing a suitable model for spin transfer 
among constituents, along with calculating higher order corrections to the various spin pro- 
cesses which can be used to test the models. There are many experiments planned at various 
locations, which will create work for the experimentalists and phenomenologists. These sug- 
gestions do not include some of the other areas of probing spin phenomena, mentioned at 
the beginning of this review. These other areas cover the spectrum of both high energy and 
medium energy nuclear physics as well. Thus, we are in an interesting period of spin physics. 
One where considerable progress has been made, only to discover that there is much more 
to be done to increase our understanding of the fundamental nature of matter through the 
use of polarization. 
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